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ABSTRACT 

In this work we study the relevance of the cosmic web and substructures on the matter and 
lensing power spectra measured from halo mock catalogues extracted from the N-body simu¬ 
lations. Since N-body simulations are computationally expensive, it is common to use faster 
methods that approximate the dark matter field as a set of halos. In this approximation, we re¬ 
place mass concentrations in N-body simulations by a spherically symmetric Navarro-Frenk- 
White halo density profile. We also consider the full mass field as the sum of two distinct 
fields: dark matter halos (M > 9 x 10 12 M 0 f. h) and particles not included into halos. Mock ha¬ 
los reproduce well the matter power spectrum, but underestimate the lensing power spectrum 
on large and small scales. For sources at z s = 1 the lensing power spectrum is underestimated 
by up to 40% at £ « 10 4 with respect to the simulated halos. The large scale effect can be al¬ 
leviated by combining the mock catalogue with the dark matter distribution outside the halos. 
In addition, to evaluate the contribution of substructures we have smeared out the intra-halo 
substructures in a N-body simulation while keeping the halo density profiles unchanged. For 
the matter power spectrum the effect of this smoothing is only of the order of 5%, but for 
lensing substructures are much more important: for { « 10 4 the internal structures contribute 
30% of the total spectrum. These findings have important implications in the way mock cat¬ 
alogues have to be created, suggesting that some approximate methods currently used for 
galaxy surveys will be inadequate for future weak lensing surveys. 

Key words: Gravitational lensing: weak - methods: analytical, numerical, statistical - cos¬ 
mology: theory, large-scale structure of Universe 


1 INTRODUCTION 


One of the most important ways to study the growth of cosmo¬ 
logical structure and infer the nature of gravity is to use gravi¬ 
tational lensing. Lensing arises from the deflection of light rays 
due to structures along the line-of-sight between the source and 
the observer. In particular, weak gravitational lensing measures the 
shear effect produced by the large scale structures in the Universe 
on the background galaxies, making it an effective tool to eval¬ 
uate the total matter i n the universe. Weak lensing observations 
are steadily improving JlU lbingei' et al. 2013; Simpsonetal . 20H 
Van Waerbeke et aP l2013l : Fu et alj2014l : Kitching et al . 2011) and 


on-going 


-beke et al.HZUliklhu et atJ20 14: Krtchmg et a 

and future surveys such as KID® DES0, Par 
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LSST0 and Euclitjf] l lLaureiis et al.l |20 1 ll : lAmendola et alJ 1 201 3l) 
will be able to probe the growth of structures and the expansion 
history to a very high precision. Such a precise determination of 
the growth of structures will provide constraints on the underly¬ 
ing cosmological model and probe whether General Relativity is 
a good description of the Universe (with cosmological constant or 
some form of dark energy) or whether modifications of gravity are 
required. 

Theoretical expectations of weak lensing on small scales 
o ften based on rayt r acing through simulations JPaceetalJ 


are otten based on raytracmg through simulatrons (Race et at. 
20071: iHilbert et al.1120091 : iFedeli et allboill: lHartlap et alj|201ll: 
Pace et al.,2011. 20151) , where light rays are shot from the observer 
to distant sources through the matter distribution of an N-body 
simulation. These simulations are useful to probe lensing statis¬ 
tics for different cosmological models and highlight differences in 
the non-linear regime that are usually not accessible via analyti- 


4 http://www.lsst.org 
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cal techniques. Many lensing studies have focused on the proper- 


convergence and shear (Munshi & Jain 

2000; Taruva et al. 

2002; 

Valageas 2000a.9; Menard et al. 20of 

Valageas & Munshi 

2004; 

Valageas et alj 2004;; Takahashi et al.l 2011), on how to model it 


how tocombinek withother cosmological probes (Songetal 


201l|: iMunshi et all 120 1 4al lbh and on tomography jMunshi et al 


2014ch . 


Much attention is also devoted to the study of 
lensing properties of a given m ass d en sity profile for 
dark matter halo s dBartelmanr] 1 996l; Qguri & Hamanal 
201l| _ |Retana-Montenegr o & Frutos-Alfarol 20 111 : 


Retana-Montenegro et al j l2012al lbl) as according to N-body 


simulations, halo density profiles can be pa r ametrised with a 
universal function (I Navarro et al.l 1 19961 . 1 19971 : iGao et al] l2008h . 
In addition, the halo density profile is jme of the main in- 
gredients fo r the halo mode l (|Schgger_&^ertschinggi 1 19911: 
ISeliakl l200d: ICoorav & Shefl] 120021: ISheth & Jairj 120031) : halo 
model application s for lensing inclu d e calculating the lensing 
power spectrum J Coorav e t all l2000l : iTakada & Hamanal 1200 3c 
ikainulainen & Marral 201 lh. the influence of substructures 
l Giocoli et ah 200S. 20101) . and the modelling of the intrinsic 
alignments dSchneider & B ridlel l20 ld) and their e ffects on the 
lensing power spectrum dCiarlariello et all |2014|) . Conversely, 
lensing can be used to constrain eleme nts of the halo model 
dCacciato et al.ll2012l : IVelander et al 120 1 4l) . 

Halo modelling also plays an important role in creating fast 
realisations of large scale structure that can play an important role 
in modelling observations of large scale structure and particularly 
their covariances. Many diff erent technique s have b een put forward 
exploi ting halo modelling. IColes & Jonesl dl99lh and ICole et ahl 
(|2005|) used the log -normal model to generate mock catalogues; 
IChuang et al .1 d2015n used an approach based on the Zel’dovich 
approximation to create mock catalogues that can accurately 
reproduce the one-p o int, tw o-po int, and thre e -point statistics; 
IScoccimarro & Shetbl d2002l) and iManera et all d2013l) used the 
2LPT formalism to create matter density fields from which halo 
catalogues were extracted (PTHalos). Ide la Torre & Peacock! 
d2013l) instead used a method based on the sampling of the mass 
function to create halos under the limit resolution of the halos; 
iTassev et alJ d20lj l2015h used a method based on Lagrangian 
Perturbation Theory (L PT) and similar methods were used b y 
dMonaco et alJl2002iJlbL iTaffoni et alJl2002l; iMonaco etal1l2013h . 

A compilation of some of th e fast methods for generating mock 
catalogues can be found in lChuang et al.1 d20l4) . Other techniques, 
based on remapping the halos from a simulation with a given cos¬ 
mology to catalogues with a different cosmolog ical background, 
have also been developed dMead & Peacock! 120131) . Finally, 
methods to generate mock catalogues on a spherical surface by 
using a real distribution were also used for the 2MASS Redshift 
Survey (2MRS) catalogue foe Domenico & Lvherisll2oT2l) . 


Halo model techniques have been tuned to reproduce galaxy 
clustering surveys, where they are exploited to calculate covari¬ 
ances of observations. It is important to understand whether they 
are similarly accurate for lensing surveys, as they can potentially 
be used to calculate not only lensing covariances, but also poten¬ 
tial cross correlations and covariances for joint lensing and galaxy 
observations. Here we attempt to answer this question. 

To do this, we create a convergence field using the halo cat¬ 
alogue extracted with a Friend-of-Friend (FoF) technique from an 


N-body simulation and compare its power spectrum with the con¬ 
vergence spectrum from the original N-body particles. We assume 
that the halos are well described by a Navarro-Frenk-White (NFW) 
density profile. In this way we can quantify the accuracy with which 
the halo catalogue reproduces the numerical simulation results, and 
on which scales we can rely on analytical methods. 

We also study the importance of substructures in the simulated 
halos. To do so, we randomly rotate the particles in a halo while 
keeping them always at the same distance from the centre. In this 
way coherent subhalos are erased while the density profile of the 
halo remains unchanged. This creates numerical halos more similar 
to the analytical ones, since the NFW profile doesn’t include the 
presence of substructures. 

Finally we explore how well fast halo catalogue techniques are 
able to reproduce full N-body simulations and the halo catalogues 
extracted from them. In particular, we focused on the PTHalos tech¬ 
nique, based on realisations of the 2LPT field ( iManera et all2013h 
and compare its results with N-body simulations. 

We begin in Sect. [2] with a brief description of the N-body 
simulations, and the implementation of the halo modelling, based 
on an NFW profile. In the following section, we compare the full 
N-body statistics to those derived from the halo catalogues, focus¬ 
ing in particular on the matter power spectrum (Sect. [3TTb and the 
lensing convergence spectrum (Sect. [T2l l and explore the impact of 
halo substructure. In Sect. [4] we investigate the effectiveness of a 
fast halo method, PTHalos, for calculating these statistics. Finally 
in Section [5] we study the effectiveness of the optimal weight on 
lensing statistics, when only mock catalogues for halos are avail¬ 
able, before drawing conclusions in Section [6] 


2 SIMULATIONS AND THE HALO MODEL 

The aim of this work is to compare results from N-body simulations 
with mock catalogues from a statistical point of view, focusing in 
particular on the matter and lensing power spectra. We begin by 
describing the N-body simulations that will be the basis of the halo 
catalogues and the methods we use to approximate them using a 
halo model approach. 


2.1 N-body simulations 

We u se the LasDama'f] N-body simulations (McBri de et al.ll2009l 
1201 lh. created with the publicly available code GADGET-II 
I Springell2005l ). In particular we use the suite of realizations called 
Oriana; the Oriana suite, an ensemble of 40 independent realiza¬ 
tions of the dark matter distribution in the Universe, studies the dark 
matter particles evolution in a box of 2.4 Gpc/h comoving filled 
with 1280 3 particles. The mass resolution is 4.573 x 10 11 M a /h 
and the corresponding softening length is 53 kpc/h comoving. Ini¬ 
tial conditions were created at z = 49 using the 2LPT formalism 
in order to take into account very early non-linearities and there¬ 
fore to be more^accurate than the usual ZeTdovich approximation 
dCrocce et al .1120061) . The in itial linear power spectrum was created 
using the CMBFAST code dSeliak & Zaldarriagal 19961) . In TableQ] 
we show the cosmological parameters adopted for the Oriana sim¬ 
ulations. 
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Cosmological Parameter 

Value 

B ary on density f2b 

0.04 

Matter density H m 

0.25 

Cosmological constant density Ha 

0.75 

Hubble parameter h 

0.7 

Spectral index n s 

1.0 

Power spectrum normalization crs 

0.8 


Table 1. Cosmological parameters adopted for the N-body simulation Ori- 
ana. 


Within the simulations, halos are identified using a paral 

"" 1985) imple- 


lel Friends-of-Friends (FoF) algorithm jDavjs_et_al 


mented within the Ntropy framework ^Gardner et al. 


2007bfal ). The 


linking length adopted is 0.2 times the mean inter-particle distance. 

In Fig. Q] we show the mass function of the Oriana simula¬ 
tion used in this work, evaluated from the FoF catalogue. We show 
the number of objects above a given mass. A halo is defined as an 
object with at least 20 particles linked together, therefore the mini¬ 
mum mass probed is M » 9 x 10 12 M e /h. As we show below, the 
sampled catalogue does not contain all of the galactic mass objects 
that might be contributing to the matter or lensing signal. 

Below, we also explore whether halo locations and masses 
derived using spherical overdensity (SO) techniques differ signifi¬ 
cantly from the Friends-of-Friends method. Flowever the halos are 
identified, all information about the internal structure of the halos 
is lost once the mock halo catalogues are constructed. One main 
goal is to quantify whether this simplification makes it impossible 
to recover the density and lensing spectra. 


2.2 The halo model approach 


To analytically model the simulated halos, we describe them as 
spherically symmetric objects characterised by a Navarr o-Frenk- 
White (NFW) density profile J Navarro et al ] |l996Lll997t) . whose 
functional form is: 


Pnfw(0 = 


Ps 

(r/r,)(l + r/r s ) 2 ’ 


(1) 


where p s = Ap c and p c is the critical density and A = 200. The 
density p s is related to the virial mass of the halo, M v ; r . The over¬ 
density is sometimes referred to the background comoving den¬ 
sity. The quantity r s is the scale radius and it is related to the virial 
radius (the radius of the halo) via the concentration parameter c, 
r s = r v „/c. In particular, the relation between p s and M vir is 


P s = 


M vi 


Anrlm{c) 


and m(c) is the dimensionless mass of the NFW halo: 


m(c) = 


I 


(i + x y 


■dx = ln(l + c) - 


1 + c 


( 2 ) 


(3) 


The only remaining unknown quantity is the concen tratio n 
parameter c. Many studies in the l i terature i Navarro et all 119971: 
Bullock et al. 200 ll : Eke et alJ 200 ll : iNeto et al.ll2007l: iDuffv et al 
20081 : Gao et al] 2008: M accio et al] 120081: Munoz-Cuartas et al 


20n|; Prada et alT 2012 ; Bhattacharva et al J 12013 ; Ludlow et al 

2013i) have been devoted to the study of the dependence of the con¬ 
centration parameter on the mass and on the redshift of the halo. 
Because low mass objects are formed earlier, when the average 



Figure 1. The number of objects above a given mass threshold M for the 
Oriana simulation at z = 0.52. 


density of the universe was higher, they tend to have higher con¬ 
centrations than higher mass objects. In this work we adopt the 
origin al concentration-mass relation established by iNavarro et al .1 
d 1997B : however, the exact recipe adopted to evaluate the concen¬ 
tration does not greatly alter our results. Differences of order of 
a few percent appear in the converge nce power spec trum on non¬ 
linear scales, with the prescription b y Ntwarroet aD dl997h giving 
more power than the iBullock et al] d2001 ) and Eke et al] 1 20011) 
recipe. 

A useful feature of the NFW density profile is that many lens¬ 
ing quantities, a nd in particular the convergence, have an anal ytical 
expression (see I B arte I m an nl 1799(3 : Iwright & Brai nerd 200(f). The 
convergence is defined as the integral, along the line of sight, of the 
density profile scaled by the critical density and weighted by a func¬ 
tion that depends on the relative distance of the lensing system com¬ 
ponents (observer, lens and source). Defining the projected matter 
density £(£?) = f Q p(0, z)dz and the critical density E c = , 


the effective convergence is therefore k(6 ) = E(0)/E crit . Here, c and 
G represent the speed of light and the gravitational constant, while 
the inverse of the lensing efficiency in the definition of the criti¬ 
cal density is given by the combination of three angular-diameter 
distances, D s , D A and D is , representing the distance between the 
observer and the source, the observer and the lens and the lens and 
the source, respectively. Finally, the vector 6 represents the angular 
position on the lens plane. 

The radial dependence of the surface mass density (justified 
by the spherical symmetry) is 


£« = 


2r s p s 
A-l) 
2r s p s 
3 

2r s p s | ] 
A-l) | 



X < 1 
X = 1 
X > 1 


(4) 


where we have defined ,v = r/r s . Note that all these equations imply 
a smooth inner profile, where substructures (density peaks inside 
the halos) do not play any role. For the comparison between the 
lensing signal of the numerical and analytical halos we will rely on 
Equation |4] and we will show its accuracy and limitations for both 
N-body (see Sect. and 2LPT simulations (see Sect. 14.41 . 


3 REPRODUCING THE N-BODY RESULTS 

In this section we study the statistical properties of the matter 
distribution in the simulated box at z = 0.52, focussing on the 
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matter power spectrum (see Sect. I3.lt and the lensing signal (see 
Sect. 13^1. 

We use a number of techniques to explore the impact of the 
substructure, and in particular how the halo catalogues perform in 
reproducing the density statistics. We directly compare the full N- 
body results to the FoF halo catalogues; to get comparable results, 
the halo catalogue must be rescaled appropriately as explained later 
so that its average density is the same as the N-body. We also di¬ 
vide the particles in the simulation into two different groups: those 
particles in halos (PIH) that, according to the FoF algorithm, build 
up the halos, and the particles outside the halos (POH) that contain 
the remaining particles. In terms of the particle numbers, it follows 
that Mot = Mpih + M>oh- This division also allows us to directly test 
the impact of treating the halos as spheres. 

To further probe the effects of halo substructure, we consider 
the impact of randomly rotating the positions of the PIH around 
the associated halo centres. We do this either by coherently rotat¬ 
ing the particle positions, such that halo substructure is preserved, 
or by randomly rotating the individual particles while preserving 
their distances to the halo centres. In the first case, the power in 
the internal structures should remain unchanged, while the cross¬ 
correlations between the substructure and the cosmic web could be 
suppressed. In the second case, the substructures are erased and the 
halos become effectively spherical, but their radial distributions are 
not fixed to a particular profile, as in the NFW halo catalogue. 

3.1 3D matter power spectrum 

The 3D matter power spectrum represents an important statistical 
tool to probe the growth of fluctuations and the underlying cos¬ 
mological model. To evaluate the matter power spectrum using the 
particle distribution in the simulated box, we divide the box vol¬ 
ume into 512 3 voxels and we assign each dark matter particle to the 
correspo nding voxel using a Nearest Grid Point (NGP) assignment 
scheme jHocknev & Eastwoodll988l) . Once the discretized volume 
is created, we Fourier transform it and multiply it by its complex 
conjugate, according to the definition of the matter power spectrum 

(5(ic)S'(k')) = (2n - $)Ps(h ■ (5) 

In the previous equation S <y> (k—k') represents the three-dimensional 
Dirac’s delta and 6 = (p - p)/p the matter overdensity. For our 
analysis, the Nyquist frequency is £ Nyq « 0.67 h/Mpc. 

A similar procedure is carried out to evaluate the mass- 
weighted matter power spectrum of the FoF halos where we assign 
the mass of each halo to the corresponding voxel determined ac¬ 
cording to the position of the halo as given in the catalogue. The 
comparison between the matter power spectrum from N-body and 
the mass-weighed spectrum is shown in Fig. [2] We also evaluate 
the shot-noise term due to the finite number of halos in the cata¬ 
logues. The halo shot noise is evaluated as P(k) = l/h, where n 
is the number density of the FoF halos. The shot-noise due to the 
finite number of particles in the simulated box is negligible. 

The two spectra agree very well on large scales where the lin¬ 
ear regime is still valid; this is obtained with an appropriate am¬ 
plitude rescaling of the halos power spectra to take into account 
the different effective matter densities of the two fields (multiply¬ 
ing the halo spectrum by the factor O m /n ha i 0 squared). Obviously 
the agreement is limited only to a restricted range of frequencies 
k < 0.1 h/Mpc, corresponding to the modes where the non-linear 
matter power spectrum starts to deviate from the linear prediction. 
For higher modes, we notice a flattening in the halo matter power 
spectrum, due to the shot noise term that dominates for high k. 



k [h/Mpc] 

Figure 2. Matter power spectrum for the N-body simulation. The red dashed 
line shows the spectrum obtained from all the particles in the simulation 
box while the blue short-dashed curve shows the rescaled mass-weighted 
spectrum of the FoF halos. 


As introduced above, we can divide the particles in the simu¬ 
lation in two different groups: the particles in the halos (PIH) and 
the particles outside the halos (POH). In Fig.[3]we show the power 
spectra for the different contributions (PIH, POH, halos and the 
cross-correlation between the PIH and the POH). For the FoF halo 
identification method described above, most of the particles are out¬ 
side the halos; the effective matter densities are £2 mPIH ~ 0.032 and 
O m ,poH ~ 0.218. As a test of consistency, we checked that sum¬ 
ming up the three spectra (particles inside and outside the halos 
and its cross-correlation) we obtain the full three-dimensional mat¬ 
ter power spectrum. 

Both the POH and the PIH show power on all scales, and are 
strongly correlated with one another. While the bias in the power 
spectrum of the halo particles is stronger, the amplitude shown is 
considerably lower since this reflects the different effective matter 
densities of each component. On linear scales, the spectra of the 
particles inside and outside the halos agree after a suitable rescal¬ 
ing of the amplitude due to the different effective matter density 
parameters. Particles in halos dominate the power seen at high k\ 
the slope in the non-linear regime for the PIH spectrum is shal¬ 
lower because the non-linear (one-halo) structures contribute to it. 
Not surprisingly, there is very good agreement between the power 
spectra evaluated from the halo catalogue and the particles building 
the halos. The agreement is excellent up to very non-linear scales 
{k x 0.2 h/Mpc) where the shot noise of the halos becomes impor¬ 
tant. 

To further evaluate the influence of the internal structure of 
halos, we rotate the particles about the halo centres as described 
above, both coherently and incoherently. The results are presented 
in Fig. [4] compared to the original PIH spectrum. When the ha¬ 
los are rotated coherently, there is very little impact on the power 
spectrum; a small decrease of power at the non-linear scales (at 
most 1%) apparently results from missing small correlations be¬ 
tween substructures in neighbouring halos. When the positions of 
individual particles are rotated inside the halos, the halo substruc¬ 
tures are wiped out and we notice a decrement in the matter power 
spectrum of order 2-3% for wavelengths k ~ 0.7 h/Mpc, the highest 
wavelength we can probe in our analysis. 

This shows two things: First, the matter power spectrum is al¬ 
most completely independent of the halo orientation; what matters 
most are the halo masses and their positions with respect to other 
halos (i.e. their correlation function). Second, the power spectrum 
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k [h/Mpc] 


Figure 3. The power spectra for the different components of the N-body 
simulation compared to the FoF halo matter power spectrum. The red 
dashed curve shows the matter power spectrum of the N-body simulation 
evaluated using all of the dark matter particles. The spectra of the particles 
in the halos (PIH) and the particles outside the halos (POH) are shown with 
violet dot-dashed and brown dot-dotted curves, respectively. The cross cor¬ 
relation power spectrum between the particles in and outside the halos is 
represented with a green dashed-dot-dotted curve. The blue short-dashed 
curve refers to the FoF halo power spectrum. 

obtained from the particles in the halos is largely independent of 
the internal substructures, just showing a decrement at the few per¬ 
cent level. This shows that most of the power in the internal halo 
distribution lies in the radial profile rather than in smaller substruc¬ 
tures. The upper and middle panel in Fig. [4] show that the power 
in halos is more impacted by the incoherent rotations than by the 
coherent rotations. This demonstrates that significant power is in ei¬ 
ther the ellipticity or the substructure of halos. Interestingly also the 
cross-correlation power spectrum is affected by the rotation of the 
particles building the halos in a similar fashion to the halo power 
spectrum. Differences are appreciable once again at small scales 
and both rotations have very similar quantitative consequences, al¬ 
beit in slightly smaller amplitude than the halo power spectrum. 
This is likely due to the induced misalignment between coherent 
substructures and the cosmic web. It is worth reminding the reader 
that on small scales, the cross-correlation between PIH and POH 
is considerably smaller than the autocorrelation of the PIH itself 
(see Fig. [3}, so the halo-cosmic web misalignment has a relatively 
small impact on the total power. While halo substructures appear to 
have an effect on the matter power spectrum limited to very small 
scales, this will not be the case any more for the lensing spectrum, 
as shown in Sect. 13.21 since this is an integrated quantity. In particu¬ 
lar, as a consequence, effects of substructures on the lensing power 
spectrum will be emphasised. 

3.2 Lensing power spectrum 

In this section we will study and compare the lensing signal from 
the full N-body simulation and from the halo catalogue as well. 

3.2.1 Creating the lensing maps 

In the following we will describe the procedure used to create the 
lensing maps taking into account the particle distribution or the ha¬ 
los in the FoF catalogue. Unfortunately, we have the matter distri¬ 
bution at a single redshift only, so we are not able to fully ray-trace 
through an evolving simulation. Since the snapshot represents the 


100000 


10000 

T 

1000 

/'" '"V v_ - 

/ i 

100 

- ! 


/ 

10 

Nbody PIH . ‘ 


Nbody, internal rotation 


Nbody, halo rotation . 


0.01 0.1 


k [h/Mpc] 



0.01 0.1 
k [h/Mpc] 


Figure 4. Upper panel: Matter power spectra with rotated halos. The red 
dashed curve shows the matter power spectrum for the N-body simulation 
using the original particle distribution in the halos. The violet dot-dashed 
curve represents a randomization of the positions of particles inside the ha¬ 
los while the brown dot-dotted curve shows the power spectrum when halos 
are rotated coherently. Middle panel: ratio between the power spectrum after 
rotations were performed and the original one. Bottom panel: ratio between 
the cross correlation power spectrum between the particles in and outside 
the halos after rotations were performed and the original one. 


simulated matter in a volume of 2.4 Gpc/h comoving, we assume 
that this same box spans the matter distribution in the universe be¬ 
tween z = 0 and z ~ 1. This is obviously not correct because we 
do not take into account the gravitational evolution of structures, 
but it is still a useful exercise to evaluate the signature of the differ¬ 
ent components under investigation, in particular how substructures 
will affect the lensing signal. 

To create the lensing maps using the particle distribution, 
we closely follow the procedure outlined in iFosalba et al.l d2008l) . 
which we refer to for more details. There are two major differ- 
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ences between our maps and the convergence maps created by 
iFosalba et alj d2008i) . Since we only have one box at a given red- 
shift, we divide the box into slices, and consider each as an inde¬ 
pendent snapshot at a different redshift and sum the contributions of 
all of them after a suitable weighting due to the lensing efficiency. 
The second difference is that we do not cover the whole sky, but 
just a portion of it; however, since lensing is a line of sight effect, 
this is not an issue. 

In the widel y _ used Bom app r oximation 

jB^elmatm^jc^eidei] 200 ll : ICoorav & Shethl [20021 : iRefregieil 
12003 ; Vale & Whitdl20 03 ). lensing distortions are integrated over 
the unperturbed photon paths and the effective convergence is 
simply defined as the weighted projected surface density: 


m = 


3 H 2 a m , 

2c 2 


£ f 


drS(f, 0) 


r(r s — r) 


( 6 ) 


where a is the scale factor, r the radial distance and r s the radial po¬ 
sition of the lensing sources. In our simulations, we assume all the 
sources are located at z s = 1, corresponding roughly to a comoving 
distance of 2.4 Gpc/h. Finally, 0 represents the angular position on 
the source and lens plane. Given the effective convergence k(0), it 
is possible to construct shear, deflection angle and lensing potential 
maps, though here we limit ourselves to the study of the effective 
convergence only. 

Since the matter distribution is discrete, we can replace the 
integ ral in Eqn. l[6]l with a summation over the lensing planes 
( IFosalba et alj2008l) . The effective convergence at a pixel (*, j ) will 
therefore be: 


x(i , j) 


3tf 0 2 n m , 0 

2 c 2 


2 S(i, j', k) 


n(r s - r k ) 

- dr k 

r s 


(7) 


where dr k represent s the plane separation (the width of the radial 
bin in lFosalba et ali I 2 OO 81 ) ). The density in each pixel, for a given 
lensing plane, p(t, j; k ), is defined as 


P(i,j',k ) 


m(i, j: k) 
dVd.j'k) 


(8) 


where dV corresponds to the volume associated with the pixel and 
m(i, j ) the mass included in the pixel (;, j) for each lensing plane 
k. Summing over all the lensing planes, we obtain the final conver¬ 
gence maps. 

We checked the numerical stability of our ray-tracing imple¬ 
mentation changing the number of lensing planes and the number 
of pixels on each plane, using either eight or sixteen lensing planes 
and a resolution of 2048 2 or 4096 2 pixels. We evaluated the lensing 
power spectrum for all the realisations and we obtained identical re¬ 
sults. This confirmed the numerical stability of our implementation. 
In the following we will therefore show results for maps obtained 
summing up 8 planes and with a resolution of 2048 2 pixels. To eval¬ 
uate the angular scale of our convergence maps, we fix a comoving 
size on the source plane of 500 Mpc/h. This corresponds, in our 
cosmology, to a field-of-view aperture of «12°, therefore the pixel 
resolution is » 21". We cut a cone through the box with side of 
0.5 Gpc/h (comoving) in order to remain in the flat-sky approxima¬ 
tion. To create more realisations of the effective convergence maps 
from our single simulation snapshot, we placed the observer in six 
different faces of the box, so that the power spectrum obtained is 
the mean of six different realisations. 

So far we have described the algorithm we implemented when 
the matter distribution is represented by the particles in the simula¬ 
tion. Something similar can be done to create the lensing maps for 
the halos of the FoF catalogue. The FoF catalogue provides for each 



Figure 5. The ratio between the lensing power spectrum with unevolving 
concentration parameter and the spectrum where the redshift evolution is 
taken into account. 


halo its mass and location, that we assume to represent the centre of 
the halo. To evaluate the convergence of the halo, we assume that 
the halos are described by a spherically symmetric NFW profile. 
Knowing the mass and the redshift of the halo, we infer the virial 
radius, the concentration parameter and the scale radius and using 
the equations described in Sect. 12.21 we evaluate the lensing signal 
as a function of the distance from the centre of the halo. Summing 
up the contribution from all the halos, we obtain the convergence 
map. As above, we verified the robustness of the implementation 
for the halos varying the number of planes and the map resolution; 
these agree very well, and to be consistent with what we did for the 
particle distribution, we use eight lens planes and a resolution of 
2048 2 pixels. 

We briefly clarify the concentration parameter we adopt for 
the halos. The concentration parameter for the NFW halo model de¬ 
pends on both the mass and the redshift of the halo; for consistency, 
when we evaluated the convergence of the halo fields, we used the 
appropriate relation discussed in Section l2/2l But in our simulation, 
the halos are all at the same redshift, making therefore the concen¬ 
tration parameter effectively dependent only on the mass. To check 
the influence of the concentration parameter on the lensing statis¬ 
tics, we constructed the convergence maps in two ways: in the first, 
we considered the concentration parameter evolving with redshift, 
while in the second case we modify the concentration accordingly 
with the mass, but assuming the same redshift for all the objects. 

We show the ratio between the two different implementations 
in Fig. [5] On large scales the two spectra are identical, while on 
small scales they differ by a few percent, with the spectrum with 
the unevolving concentration being slightly higher. Larger differ¬ 
ences appear at intermediate frequencies, where the spectrum with 
unevolving concentration parameter is roughly 8% higher than the 
spectrum where the concentration parameter is a function of mass 
and redshift. In the following we will assume that the concentra¬ 
tion varies with mass and redshift and will discuss the impact of 
this assumption. 

Another comment is necessary at this point. Comparing our 
results to better ray-tracing simulations, we need to emphasize that 
there will be additional contributions. Not only the concentration 
relation will evolve in time, but also the mass function will play 
an important role in changing the number of halos as a function of 
redshift. 
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Figure 6. Effective convergence power spectrum for the total matter distri¬ 
bution (red dashed line) and the FoF halo catalogue (blue dotted line). 

3.2.2 Comparing the convergence spectra 

In Fig.[6]we compare the N-body tensing power spectrum, obtained 
using the particles in the simulated box, with the spectrum obtained 
using the halos only. The lensing power spectrum is significantly 
different from the total spectrum; unlike the matter power spec¬ 
trum, they cannot be matched with a simple rescaling to account 
for the effective matter density. The two spectra have very different 
shapes. The lensing spectrum of the halos has a lower amplitude 
(since f2 m>h < fi m ) and differs both at low and high £s, while it is 
comparable around { ~ 2x 10 3 . This suggests the importance of 
the cosmic web - structures not in halos - for the lensing signal, be¬ 
ing them the primary responsible for the differences of shape and 
amplitude on large scales for the spectrum obtained using only the 
halos. 

As above, it is useful to break down the contributions of the 
particles between those making up the halos (PIH), and those par¬ 
ticles outside the halos (POH). In Fig. [7] we show the contribution 
to the total lensing spectrum of PIH, the POH and of their cross¬ 
correlation. Again, we verified that summing of the spectra of the 
different components recovers the power spectrum of the full field. 

The lensing spectrum inferred from the NFW halo catalogue 
agrees very well with that of the PIH on large scales, up to t » 200, 
but there is an excess of power around C ~ 1000 of about 10%. 
However, on small scales we observe a lack of power inferred from 
the halo catalogue; we find the power dips about 20% at t ~ 10 4 , 
dropping even lower on smaller scales. The origin of this complex 
behaviour is unclear; it must be remembered that the lensing maps 
are a projection of the density maps, and the effects of substruc¬ 
ture missing in the halo treatment (as well as associated shot noise) 
could be spread over a range of multipoles. 

Different concentration recipes give very similar results on 
small scales, but they can impact the excess of power seen on in¬ 
termediate scales. The peak seen in the power ratio can change as 
much as 5% and shifts towards smaller frequencies when chang¬ 
ing the concentration relation. When considering FoF halos, the 
amplitude and the position of the excess is related to the assumed 
concentration relation, but its origin is more complex and related to 
the asphericity of the FoF halos, as we will explain later. 

In Fig.[8]we compare the lensing signal of the full simulation 
to that combining the NFW halos with the matter distribution of 
particles outside the halos (upper panel), showing their ratio in the 
lower panel. These are largely in agreement, but are unable to re¬ 
produce the signal exactly, due to the lack of power on large t s in 
the NFW halos and the mismatch on intermediate scales (see lower 



Figure 7. Upper panel: convergence power spectrum for the different com¬ 
ponents in the N-body simulation and the NFW halo catalogue. Lower 
panel: ratio between the power spectra of the halos and of the PIH. Colours 
are as in Fig. [7] 


panel in Fig. [Til. Due to the importance of the POH and its cross 
correlation with the halos, the lack of total power is limited to 10% 
and changes little over a decade in scale. 

We also checked how these results are affected by the defi¬ 
nition of a halo itself, by analysing maps created with catalogues 
whose halos had a different particle number threshold. All struc¬ 
tures having fewer particles were excluded from the halos and the 
particles building them were absorbed into the cosmic web. Thresh¬ 
olds of 20, 100, 200, 1000 and 2000 particles, corresponding to a 
mass threshold of 9 x 10 11 , 4.6 x 10 12 , 9 x 10 12 , 4.6 x 10 13 and 
9 x 10 L ’ M 0 //i, were considered. The differences between the com¬ 
bined POH plus mock catalogues and the full particle distribution 
for the halos become progressively smaller when fewer halos are 
used (and the fraction of mass in the POH increases). 

One might wonder whether the differences using the analyti¬ 
cal halos could be related to how the halos were identified in the 
simulations. The NFW profile assumes that halos are spherically 
symmetric objects with a given over-density (in our case 200 times 
the critical density). Obviously this is not necessarily true in the 
FoF halos that they replace. To explore whether the spherical sym¬ 
metry of the halos is an issue, we apply a spherical over-density 
algorithm on top of the FoF halo catalogue to identify the particles 
in the halos. The halos identified in this way are by definition more 
spherical; we also require the number of particles must be at least 
20 and the density enclosed in the sphere to be p = 200p c . These 
spherical over-density (SO) halos will be smaller (and therefore less 
massive) and in a lower number with respect to the FoF halos. To 
create the SO-halo catalogue from the FoF-halo catalogue, we con- 
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Figure 8. Upper panel: convergence power spectrum from the full particle 
distribution (red dashed line) and from the combination of the cosmic web 
(POH) and the NFW halos (grey solid line). Lower panel: the ratio between 
the spectra. 


sider a sphere, centred on the halo centre with radius equal to the 
maximum distance of the particles from the centre of the halo; we 
then shrink this radius until the internal average density is 200 times 
the critical density. If the final number of particles is less than 20 
(the minimum number of particles to have an halo) we discard that 
halo and any discarded particles instead contribute to the cosmic 
web. 

In Fig. [9] we show the comparison between the lensing signal 
from the particles building the numerical SO halos and the corre¬ 
sponding analytical NFW halos. The agreement on large scales is 
once again achieved, as expected. Up to £ ~ 400 the two spectra 
are practically identical while with the increase of the multipole we 
notice the usual lack of power, while we find the power decreases 
rather than increases at intermediate scales. Comparing Fig.[9]with 
Fig.[8]we notice that the previously discussed peak disappears. This 
tells us that the FoF halos are responsible for excess of power at in¬ 
termediate scales. When we compare a spherical halo with a FoF 
halo, we observe that the latter will have particles outside the mod¬ 
elled spherical area, creating overdense regions that by construction 
do not exist in the case of the SO halos. These regions are respon¬ 
sible for the excess of power at intermediate scales. 


3.2.3 Rotating the halo particles 

We next explore the impact of rotating halo particles, either ran¬ 
domly or coherently, on the lensing convergence spectrum. Recall 
that for the power spectrum, these rotations had relatively little im¬ 




Figure 9. Upper panel: convergence spectrum from the particles building 
the SO halos and the corresponding NFW ones. Lower panel: ratio between 
the spectra. 


pact; however, by suppressing halo substructures, the incoherent 
rotations did suppress the large k spectrum (Fig.0. 

The impact of the rotations for the lensing spectrum is shown 
in Fig. [TO] There is no appreciable difference on large scales; both 
kinds of rotations remain in very good agreement with the origi¬ 
nal PIH spectrum. However, on small angular scales, where one- 
halo effects become important (above l » 1000), there is a notable 
lack of power for the incoherent rotations that becomes as large 
as x; 30% before the shot noise contribution becomes dominant at 
l ss 10,000. This indicates that the halo substructures play a larger 
role in the determination of the lensing maps than they do for the 
density maps, due to closer substructures being projected into larger 
angular scales. It may be possible to improve the halo model results 
by adding substructure t o the assumed profi l es to reflect the mass 
and profile distributions dGiocoli et al.ll200i3 . l20 id) of the satellite 
sub-halos. Alternatively, one could create a library of halos includ¬ 
ing substructures and assign them randomly to the mock halos. 

It is interesting to compare these results to the case where the 
particles in a halo are rotated coherently, preserving the halo sub¬ 
structures. In this case (also shown in Fig. 1 1 0b. the power is largely 
preserved compared to the original simulations (as also seen in the 
matter power spectrum); this indicates that it is indeed the halo sub¬ 
structure, rather than the halo profile, that makes the differences 
seen for the incoherent rotations. Mixed in with the effect of sub¬ 
structure will also be the ellipsoidal distortions seen in realistic ha¬ 
los; while spherically symmetric objects have no preferred direc¬ 
tion, the effective convergence of the halos (ellipsoidal rather than 
spherical) will change significantly if the projection is made along 
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Figure 10. Upper panel: Comparison between the effective convergence 
power spectra for the original simulation (red dashed curve) limited to the 
particles in the halos, and the same particle distribution but with a rotation 
about the centre of each halo (violet dot-dashed curve) or with randomiza¬ 
tion of the particles’ positions (brown dot-dotted line). Lower panel: ratio 
between the different rotated matter distribution spectra and the original 
one. The red dashed and blue dot-dotted lines represent the spectra of the 
NFW halos, as a ratio to the FoF or SO simulated halos, respectively. 


the longest or shortest axis. Thus the ellipsoidal structure of halos 
could add additional variance in the lensing power spectrum. 

From the analysis of Fig. m we can draw another impor¬ 
tant conclusion. By comparing the curves showing the ratio be¬ 
tween the lensing power spectrum based on the FoF halo catalogue 
(red dashed curve) and on the SO halo catalogue (blue dot-dotted 
curve), we notice that the peak at { * 1000 disappears, clearly 
showing that its presence can be related to the non-spherical halos 
typical of FoF catalogues as well as bridged halos (unbounded par¬ 
ticles included by the FoF algorithm). We want to stress out that the 
peak at t » 1000 is the effect of using FoF halo catalogues instead 
of SO halo catalogues on lensing maps. 

Finally, we investigate how the substructure contribution de¬ 
pends on the masses of halos we consider. In Fig. |TT] we investigate 
how the effect of substructures is affected by the different definition 
of minimum halo mass considered, plotting the ratio of the effective 
convergence power spectrum of halos whose internal substructures 
are washed out to the PIH power spectrum for different halo mass 
thresholds. Massive halos are rarer but are rich in substructures and 
contribute more to the total halo mass. Thus, it is not surprising that 
erasing substructures in massive halos leads to a bigger suppression 
than when all the halos are included. Changing the threshold by a 
factor of fifty doubles the fraction of suppressed power and shifts 
it to larger angular scales. Given the rarity of very massive objects, 



Figure 11. Ratio between the power spectra of the halos with internal sub¬ 
structures washed out and the original convergence power spectrum of the 
numerical halos for different particle number (mass) threshold. Different 
colour and line styles correspond to different particle thresholds. The red 
dashed line corresponds to halos with at least 20 particles; the blue dotted 
(violet dot-dashed) line to halos with at least 100 (200) particles; the brown 
dot-dotted (grey dashed dot-dotted) line to halos with at least 1000 (2000) 
particles. A threshold of 20, 100, 200, 1000 and 2000 particles correspond¬ 
ing to a mass threshold of 9 X 10 11 , 4.6 X 10 12 , 9 X 10 12 , 4.6 X 10 13 and 
9 X 10 13 M B /h, respectively. 


we do not see appreciable differences with a threshold of 1000 or 
2000 particles. 

Our inability to reproduce the full_simulation results should be 
compared with the findings of IPetkova et al.l d201 3l) . In that work, 
the authors present an ext ension of the raytracing code GLAMER 
jMetcalf & Petkov3l201 3|> to the case of multiple deflection planes 
and compare their findings with several existing works. In partic¬ 
ular they assemble spherically symmetric NFW halos drawn from 
the Sheth & Tormen mass function and compare the probability 
distribution functions for deflection angle, shear and effective con¬ 
vergence to N-body results from Pac ket all l201ll) . finding good 
agreement. However, in IPace et alj 1 21)1 5t) it was shown that for 
high values of the effective convergence, the raytracing simula¬ 
tions underestimated analytical or semi-analytical solutions. Thus 
the agreement may have been due to the suppressed power observed 
here for NFW halos compensating for the loss of resolution in the 
procedure adopted to build the lensing planes (for a more complete 
discussion see [Pace et al .203). Note that, because we directly use 
the particle distribution to evaluate the effective convergence rather 
than evaluating it as the Laplacian of the lensing potential, we ex- 
pect the resolution issues her e to be less severe than in IPace et al.l 
00071) . [Pace et al.l j201ll) and lPace et alJ d2015h . 

To conclude, our results show that both the cosmic web and 
the halo substructures are necessary to fully reproduce the lensing 
statistics, suggesting that halo based approaches may be seriously 
flawed for lensing studies. 


4 2LPT METHODS 

In this section we focus on the comparison between N-body simu¬ 
lations and simulation and halo catalogues based on the Lagrangian 
Perturbation Theory (LPT). 
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4.1 Lagrangian perturbation theory 


4.2 2LPT halo catalogues 


Lagrangian perturbation theory (LPT) is a very useful mathemat¬ 
ical approach used to describe perturbations in the matter density 
field and to create initial conditions for cosmological N-body sim¬ 
ulations lBuchert||l989|; Moutarde et al.|l99l|; lBouchet et all 

"iBucht 


Buchert 1992; Buchert & Ehlers 1993; Buchert 1994; H ivon e t alJ 


dy sim- 

H 19921; 


19951 ; Catelan|[l995l) . In the LPT formalism, the current Eulerian 

position fofa given mass element is related to its initial Lagrangian 
position 4 via a displacement vector field 'i(q) according to the re¬ 
lation 


2=4+ ^(q) ■ 


(9) 


Similarly to what it is usually done for the Eulerian perturba¬ 
tion theory (EPT), one can make a power series expansion of the 
displacement field 'f as a function of the linear over-density field 

S L : 


Y = vpd> + >p <2) + >P (3) + 'P <4) + ... , (10) 

and 'F (n) is related to the /7th order expansion of the linear mat¬ 
ter density field. LPT up t o third order has been d i scussed i n sev¬ 
eral works in literature (Buchert & Ehler j 1 19931 ; iBuchertl 1 1994 
iHivon et all 1 1 9951 : ICatelanl 1 1 995h and recently cal culations at the 
fourth order were carried out by iTatekawal <20131) . Here, we will 
limit ourselves to the second-order treatment. Translating from Eu¬ 
lerian to Lagrangian coordinates, we can obtain the LPT equations 
at each order by substituting Eqn. ( 1101 into Eqn. <(9} and the equa¬ 
tion of motion for particles 


tPl dj? 1 

-t+2H— =--V0. 

dt 2 dt a 2 


( 11 ) 


where a is the scale factor, H the Hubble expansion rate and (j> the 
gravitational potential. 

At first orde r we obtain the well-known ZeTdovich approxi¬ 
mation feel’ Dovichll 1 970l) : 

= ( 12 ) 


In Eqn. 412b . D\(a) represents the linear growth factor and 6(4) 
is the Gaussian density field imposed by the initial conditions. 
Second-order LPT (2LPT) corrects the ZeTdovich approximation 
by taking into account gravitational tidal forces and the second or¬ 
der displacement V P ,2) to find 


- ■ = \d 2 (o) £ . 'I'. - T„ 'I',,) , (13) 


i*i 


where D 2 is the second order growth f actor and for a wide range of 
£2 is proportional to D\ (Bouchet et al.il 19951) . 

At first and second order, it is convenient to define two dif¬ 
ferent gravitational potentials 0 (l) and 0 (2) since the solutions are 
curl-free. The displacement fields are related to the gravitational 
potentials via the relation 

, (14) 

with n = 1,2. Using the relations in Eqn. [T4] the solution up to 
second order is 

m = (is) 

* = + , (16) 

where v = is the peculiar velocity and t the cosmic time. Finally, 
f is the logarithmic derivative of the growth factor, f = . 


With the formalism outlined in Sect. 14.11 it is possible to gener¬ 
ate the 2LPT displacement and evolve particles according to it. 
From these, halos can be identified and rescaled to provide fast 
h alo catalogues (PTHa l os); we briefly descri b e these here and refer 
to lManera et al.1 d2013l) : IChuang et all d2QI5t) : lManera et alj d2015h 
for more details. 

The algorithm used to create the second order displacement 
V P <2) is described in IScoccimarrol d 19981) : it takes advantage o f 
the fast Fourier transform ( FFT) technique. IColes et ail d 1 9931 ); 
IScoccimarro & Shethl d2002h used a cut-off in the linear matter 
power spectrum and found the halos using merger trees. For these 
simulations halos are identified with a standard friends-of-friends 
(FoF) technique dDavis et alj 19851) with a modified linking length. 
The theoretical derivation for the evaluation of the modified value 
for the 2LPT linking length is based on the spherical collapse model 
for the 2LPT and we refer to lManera et al .1 42013|) for an exhaustive 
explanation. The linking length used for the 2LPT simulations is 
calibrated on the LasDamas N-body simulations and its value is 
b = 0.38, significantly larger than the standard b = 0.2 commonly 
used in N-body simulations. 

The FoF algorithm identifies halo positions and masses in the 
2LPT realisations. However, the 2LPT masses need to be calibrated 
by abundance matching to match a fiducial mass function. In par¬ 
ticular, the lTinker et all 120081) mass function was used to calibrate 
the 2LPT halo masses. 

As for the LasDamas simulations used previously, 2LPT halos 
are identified at z = 0.52. The cosmology adopted is, as before, a 
flat ACDM model with £2 = 0.25, £2 A = 0.75, £2 b = 0.04, h = 0.7, 
power spectrum normalization erg = 0.8 and spectral index n s = 1. 
The simulated box is 2.4 Gpc/h comoving on a side with 1280 3 dark 
matter particles. In the following we will consider objects with at 
least 20 particles. 


4.3 3D matter power spectrum 

We first study the matter power spectrum following the same pro- 
cedur e outlined in Sect,l3.1| As fo r the full simulations, we use the 
NGP dHocknev & Eastwood! 1988i) algorithm to assign dark matter 
particles to the grid and we use 512 points for the Fourier transform. 

In the upper panel of Fig.[l2]we first compare the full N-body 
spectrum to that inferred from the 2LPT particles, showing their 
ratio. As expected, the 2LPT spectrum agrees very well with the 
N-body spectrum on large scales, but the spectra differ consider¬ 
ably on small scales where fully non-linear effects are absent in the 
2LPT treatment. The power deficit is considerable on the smallest 
modes of the simulation, dropping to 20% of the full simulation. 
As we will see in the next section (see Sect. l4.4| i. this inaccuracy 
will have tremendous effects on the lensing power spectrum. 

It is necessary to find alternative methods to enhance the small 
scale power in the 2LPT simulations. In the halo method, the under¬ 
developed 2LPT halos are replaced by more realistic halo masses. 
As pointed out bv lManera et al] d2013l) . to correctly reproduce the 
fiducial mass function, it is necessary to calibrate the mass of the 
FoF halos. We therefore compare the halo power spectrum with the 
total matter power spectrum in the bottom panel of Fig. \n\ where 
we refer to the caption for the different line styles and colours. As 
done for the FoF halos in the N-body simulation, we evaluate the 
halo power spectrum weighting the halos either by their (corrected) 
mass, or simply giving all equal weight. 

The halo power spectra reproduce, to large extent, the matter 
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Figure 12. Upper panel: Ratio between the 2LPT and the N-body matter 
power spectra (blue short-dashed curve). Lower panel: Comparison of the 
matter and halo power spectrum for the 2LPT simulation. The red dashed 
and blue short-dashed curves represent the N-body and the 2LPT matter 
power spectra, respectively, while the brown dot-dotted curve represents 
the halo power spectrum with mass weight. 


Figure 13. Upper panel: Effective convergence for the N-Body simula¬ 
tion (red dashed curve) and the 2LPT simulation (blue short-dashed curve). 
Solid line shows the shot-noise power spectrum. Lower panel: ratio between 
the particle-based 2LPT simulation and the N-body simulation. 


power spectrum on large scales. We checked that rescaling the halo 
power spectrum by an appropriate factor, well reproduces the ex¬ 
pected matter power spectrum. The halo power spectrum with mass 
weighting is indeed very close to the total one due to a rescaling in 
the matter density parameter. 

4.4 Lensing power spectrum 

We next analyse the lensing power spectrum in the 2LPT simula¬ 
tions, again beginning with the 2LPT particle distribution and then 
seeing how things change using FoF halos derived from the simu¬ 
lations. 

In Fig. [I3]we show the lensing power spectrum using all the 
particles in the simulations for sources at z s = 1 for both the N- 
body and the 2LPT simulation. The missing power seen in the mat¬ 
ter power spectrum (due to the lack of non-linear evolution in the 
2LPT simulations) is even more prevalent in the effective conver¬ 
gence power spectrum. For very low t the agreement is still very 
good, but the power deficit is significant on intermediate scales, be¬ 
ginning as low as C ss 100. By £ « 1000, the power drops to 20% of 
the full simulation; this is as bad as the deficit seen in the density 
spectrum only at the highest k modes. Beyond C ~ 1000. the lensing 
spectrum flattens out at the expected shot nose level. It is perhaps 
not surprising that the missing non-linear power of the 2LPT simu¬ 
lations is disastrous for predictions of the lensing statistics. 

Better predictions might be derived using the associated NFW 
halo catalogues, where the halo masses have been rescaled to match 


the mass distribution observed in full N-body simulations; boost¬ 
ing the halo masses should make the deficit shown in Fie.ll3lless 
severe. In Fig. [T4] we show the spectrum of the PTHalos plus the 
POH in the cosmic web. We indeed see that the power decrement of 
the 2LPT particles has disappeared, but the mass rescaling makes 
the signal too high over a range of frequencies. While on large 
scales the agreement is excellent, on scales with t > 100 we instead 
note an excess of power with respect to the total power measured 
from the N-body simulation. This is due to the mass correction we 
applied to the halos before evaluating their lensing signal, so that 
the mass of the rescaled halos added is greater than the PIH initially 
removed. This shows that PTHalos are clustered differently. 

In Fig.[l5]we compare the effective convergence power spec¬ 
trum of the halo catalogues for the N-body and the 2LPT simula¬ 
tions, taking into account the mass correction for the 2LPT halos. 
While the mass correction is able to populate the halosywith galax¬ 
ies following the HOD prescription (Manera et alj|2013i) . the ef¬ 
fective convergence for halos from the 2LPT catalogue is initially 
lower than the spectrum of the halos from the N-body catalogue 
on all scales. However, the total mass in the N-body halos is still 
greater (» 30%) than that in the 2LPT halos even after the masses 
are individually rescaled. When the effective convergence power 
spectrum is rescaled by the matter density parameter squared, there 
is good agreement with the N-body lensing spectrum, as shown 
in the figure. However, the spectra do not match over the whole 
range of frequencies for the halo lensing power spectrum; there re¬ 
mains a difference in shape, with the N-body having more power 
at t ~ 2000. This might reflect the fact that the mass remapping is 
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Figure 14. Convergence power spectrum from the N-body full particle dis¬ 
tribution (red dashed line) and from the combination of particles building 
the cosmic web and spherical symmetric halos (grey solid line) of the 2LPT 
simulation. 



Figure 15. Effective convergence for the N-body halos (red dashed curve) 
and the 2LPT PTHalos (blue short-dashed curve). Halos are described by a 
spherical symmetric NFW density profile. 

not perfect, leading to a somewhat different mass weighting of the 
clustering of the PTHalos than in the N-body simulation. 


5 INCLUSION OF A BIAS WEIGHT 

Given the limitations of simply using the mock-catalogues when 
trying to reproduce the results from the N-body simulations, in par¬ 
ticular those for weak lensing, we next investigate whether these 
issues can be alleviated by using an optimal weight to reconstruct 
the underlying dark matter distribution from the halo c atalogue. To 
do so, we closely follow the work by ICai et al.l ( |201 ll) and we re¬ 
fer to them for further details. The idea is to see whether lensing 
studies can benefit from a mass-dependent weighting used to re¬ 
construct the three dimensional matter field starting from a mock 
catalogue. 

The halo distribution represents a biased sample of the under¬ 
lying total matter distribution <5 m and the goal is to reconstruct the 
density only using the information we can extract from the mock 
catalogue. We divide our mock catalogue into mass bins and indi¬ 
cate with <5, the fluctuation in each bin. The estimator <5 m for the 
matter density field is then 

S m = ^ w t S, , (17) 




k [h/Mpc] 

Figure 16. Optimal weight for the five different mass bins as a function 
of the wavelength k. From top to bottom we show the optimal weight for 
bin 1 (red dashed line), 2 (blue dotted line), 3 (orange dot-dashed curve), 
4 (purple dot-short-dashed curve) and 5 (brown dot-dotted curve). Upper 
(lower) panel shows the weights for the N-body (2LPT) simulations. 

where w, represents the weight we use for each mass bin to recon¬ 
struct the dark matter field starting from the mock catalogue. The 
power spectrum for the mass fluctuation P, the bias vector b and 


the covariance matrix C are defined as 

<<0 = P, (18) 

<<5nA> = b,P, (19) 

(8fij) = Cij . ( 20 ) 

Our weights are determined by minimising the stochasticity 
between the true and reconstructed density fields. Defining the 
stochasticity E as 

E 2 = l-2b T w+ w T Cw/P , (21) 

the minimum in the stochasticity is reached when the weight is 

WV = (C/PT'b , (22) 

this implies the optimal stochasticity is 

E 2 0ft = l-b T (C/P)'b . (23) 


The halos are divided into five (logarithmically) equally- 
spaced mass bins, from the smallest masses (halos with 20 parti¬ 
cles) in bin 1, to bin 5 which includes halos up to 70% of the highest 
mass found in the simulation. This most massive bin has only nine 
halos. To evaluate the weights, we calculated the power spectrum 
for each bin and the cross correlation between the bins and the total 
matter field. 

In Fig.[l6]we show the optimal weights to be applied to each 























Cosmic web and substructure for spectra 13 



k [h/Mpc] 



k [h/Mpc] 

Figure 17. Comparison between the total mass power spectrum (red dashed 
line) with the reconstructed one using the optimal weights evaluated in 
Eq. 1 22 1 (blue short-dashed line). Upper (lower) panel shows results for the 
N-body (2LPT) simulations. 


vergence k(9), for sources at a fixed redshift, is written as 



, , hiw'if. k(w ~ W) S[f K (w'6), w'] 

aw - 

Jk(w ) a 


(24) 

where a is the scale factor, w the comoving distance and f K a func¬ 
tion that takes into account the 3D curvature of the space-time. In 
Fourier space, 6 can be viewed as the weighted sum of the contri¬ 
butions of all the different mass bins (see Eq. 117b . Therefore the 
effective convergence power spectrum is the integration along the 
line of sight of the reconstructed matter power spectrum: 


M)in 


M)in 


P&) = 2 w2 d k ) p m + 2 ^ wMwjOoP.jdc) 


i=l,j>i 


(25) 


where P,(k) represents the matter power spectrum of the i-th bin 
and Pjj(k ) the cross-correlation between the different bins. Using 
these relations, we numerically evaluate the lensing power spec¬ 
trum from the halo catalogue. The matter power spectrum is binned 
only for a given set of wavelengths, therefore when the evaluation 
of the matter power spectrum has to be done outside this range, we 
extrapolate the values of the spectrum to both smaller and larger 
scales. 

We show our results in Fig. |T8] As expected the reconstruc¬ 
tion is not perfect and is limited to a limited range of frequencies. 
On large scales the reconstructed lensing power spectrum agrees 
well with the original one, but at higher wave numbers the recon¬ 
structed spectrum has a lower amplitude: this reflects the lack of 
power present in the reconstructed three-dimensional matter power 
spectrum. Interestingly, the reconstruction for the 2LPT case works 
better than for the N-body simulation. However, what is recon¬ 
structed is the lensing prediction for the 2LPT particles, and so 
greatly under-produces power compared to the N-body lensing pre¬ 
diction, as seen in the bottom panel of Fig.[l8l 


mass bin. Due to the limited numbers of objects, the weights are 
very noisy on large scales, particularly for the highest mass halos. 
It is easy to understand why the lowest mass bin has a higher nor¬ 
malisation: lower mass halos are less biased and better reflect the 
matter density field not in halos, hence they are up-weighted with 
respect to the other halos in the density reconstruction. This can 
be seen in both the N-body and the 2LPT simulations. On small 
scales, the amplitude drops a factor of 2.5 going from the lowest 
mass bin to the highest mass bin in the N-body simulation, while 
for the 2LPT simulation differences in the weights are more pro¬ 
nounced. 

In Fig. [HI we show the reconstructed three dimensional mat¬ 
ter power spectrum using the mass weighting. The reconstruction 
on large scales is relatively noisy, though the reconstructed mat¬ 
ter power spectrum follows the expected trend. Good agreement is 
achieved for the mildly non-linear regime, with differences of the 
order of 10% for k ~ 0.15 h/Mpc. For smaller scales, deeply in the 
non-linear regime, we observe a lack of power that becomes more 
pronounced for higher wave numbers. The effect of the shot noise 
is suppressed with the weighting, making the spectrum more reli¬ 
able over a larger range of frequencies. In addition, now the spec¬ 
tra agree over a large range of frequencies without any additional 
rescaling due to the different effective matter density parameter, as 
the weights compensate for this. 

In terms of the matter density perturbations, the effective con- 


6 CONCLUSIONS 

Our focus has been to understand how well approximations to large 
scale structure, in particular those based on the halo model, can re¬ 
produce the statistics characterising large scale structure observa¬ 
tions. Particularly challenging in this regard are lensing statistics, 
where projection effects mean that large and small scale features 
can jointly determine the observations on a given angular scale. 
Lensing studies are often based on ray-tracing through numerical 
simulations, but these methods are computationally expensive and 
it is essential to understand which fast methods can make reliable 
lensing predictions. 

In this work we have studied whether mock halo catalogues, 
based either on full simulations or fast 2LPT realisations, could be 
of value for pushing lensing analyses into the non-linear regime. 
Such halo model approaches, if based on simple templates such as 
the NFW profile, can take advantage of analytical methods to cal¬ 
culate the lensing effect and so be much faster. By comparing fully 
non-linear simulations to mock catalogue methods over a huge cos¬ 
mic volume, we have quantified where the halo model approaches 
appear to fail. 

We have considered a gradual progression from the full N- 
body simulations to the halo catalogues. The first step beyond di¬ 
rectly using the N-body simulation was to coherently rotate the ha¬ 
los; and we found that the halo orientations did not have a signif¬ 
icant impact on either the matter power spectrum or the lensing 
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Figure 18. Comparison between the tensing power spectrum obtained from 
the full the matter distribution (red dashed line) and the lensing power spec¬ 
trum evaluated from the reconstructed density field (blue short-dashed line). 
Upper (lower) panel shows results for the N-body (2LPT) simulations. 


convergence spectrum. However, when the halo particles were ro¬ 
tated incoherently, effectively smearing out the internal halo sub¬ 
structures, the matter power spectrum was reduced by 5% on the 
smallest scales we probed. The lensing convergence was impacted 
more dramatically, as this blurring of halos resulted in a signifi¬ 
cant suppression of power even on intermediate angular scales (see 
Fig.®. 

Directly replacing the halo particles with NFW halos had a 
similar effect, though in detail it depends on the assumed prescrip¬ 
tion for the concentration evolution. We also find that the effect 
of ignoring substructure is stronger in bigger halos. Finally, when 
we consider the halo catalogues alone, we can still reproduce the 
density power spectrum on large scales with a suitable rescaling, at 
least up to scales k » 0.1 h/Mpc, where shot noise in the catalogue 
begins to dominate. But this cannot be done for the convergence 
spectrum; on large scales, the cosmic web of particles not included 
in halos plays an essential role in determining not only the ampli¬ 
tude but also the shape of the lensing spectrum. 

One fast means of creating halo catalogues is through a 2LPT 
approach. 2LPT realisations entirely miss fully non-linear struc¬ 
tures, and on their own would be inadequate for pushing predictions 
into the non-linear regime. Methods such as PTHalos improve the 
non-linear properties by identifying embryonic halos in the 2LPT 
simulations and replacing them with more realistic halo masses 
based on abundance matching to full simulations. These can im¬ 
prove agreement in the halo statistics, including the lensing power 
spectrum. 

We also considered an optimal reconstruction of the matter 


density field by reweighing halos according to their mass, follow¬ 
ing the work bv lCai et al.1 11201T). In this case the reconstructed mat¬ 
ter power spectrum agrees reasonably well with the original matter 
power spectrum obtained from the particle distribution and no scal¬ 
ing is necessary to compensate for the different matter density pa¬ 
rameters. The method worked better on large scales in reproducing 
the shape of the lensing convergence spectrum, compared with the 
simple halo approach. However, the lack of power on small scales 
remains an outstanding issue. 

We demonstrate that both cosmic web and substructures are 
important in shaping the total lensing signal. While the cosmic weh 
is responsible for the overall shape and contribution on large scales, 
substructures are responsible for a quite substantial fraction of the 
power on small scales. 

To improve the performance of halo methods for lensing on 
non-linear scales, at a minimum we must include some additional 
substructure beyond the assumption that halos have simple spher¬ 
ically symmetric profiles. An obvious possibility to explore is the 
impact of the halo ellipticity; it would be interesting to understand 
whether extracting and modelling the halo ellipticities could help 
address the power deficits we have seen in the lensing spectrum. 
Another possibility is to include substructur es, suc h as s atellite 
galax ies, following, for example, the work of iGiocoli et alj J2008L 
l20ld) in characterising the number and masses of satellite halos. 
Alternatively, one could create a library of realistic halos (or their 
associated convergence maps) that could be pasted randomly onto 
halo catalogues. 
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